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A  FAST  CONVERGING  CASCADED  ADAPTIVE  CANCELLER 


1.  Introduction 

The  use  of  adaptive  linear  techniques  to  solve  signal  processing  problems  is  needed  particularly  when 
the  interference  environment  external  to  the  signal  processor  (such  as  for  a  radar  or  communication 
system)  is  not  known  a  priori.  Due  to  this  lack  of  knowledge  of  an  external  environment,  adaptive 
techniques  require  a  certain  amount  of  data  to  cancel  the  external  inteiference.  The  amount  of  data  (the 
number  of  independent  samples  per  input  sensor)  required  so  that  the  performance  of  the  adaptive 
processor  is  close  (nominally  within  3  dB)  to  the  optimum  is  called  the  convergence  measure  of 
effectiveness  (MOE)  of  the  processor.  The  minimization  of  the  convergence  time  is  important  since  in 
many  environments  the  external  interference  changes  rapidly  with  time. 

With  assumption  that  the  input  interference  is  Gaussian,  the  classical  adaptive  linear  processor  with 
N  inputs  is  based  on  forming  the  sample  covariance  matrix  via  maximum  likelihood  (ML)  estimation.  The 
linear  weight  vector  is  found  by  multiplying  the  inverse  (assuming  it  exists)  of  the  sample  covariance 
matrix  by  a  desired  steering  vector.  The  convergence  MOE  of  this  technique  which  is  called  the  sample 
matrix  inversion  (SMI)  algorithm  [1]  is  (with  some  assumptions)  independent  of  the  external  noise 
environment  and  is  approximately  twice  the  number  of  independent  sensor  inputs.  However,  if  fewer 
samples  are  available,  performance  will  degrade.  For  example,  if  there  were  only  one  narrowband  (NB) 
jammer  present  and  N  —  50  sensors,  the  SMI  requires  roughly  K.  2A  100  samples  per  input 
channel.  Intuitively  only  a  few  samples  should  be  required  since  only  one  NB  Jammer  is  present. 

Due  to  the  failure  of  the  existence  of  the  ML  solution  ior  K  <  N  samples  (the  sample  covariance 
matrix  is  singular),  there  have  been  several  techniques  proposed  to  improve  convergence.  Techniques 
such  as  loaded  sample  matrix  inversion  (LSMI)  [2],  the  adaptive-adaptive  technique  [3],  and  subspace 
techniques  [4,5,6]  have  been  proposed.  However,  a  number  of  these  techniques  are  based  on  heuristic 
constructions.  Most  subspace  techniques  such  as  in  [4,5,6]  are  based  on  utilizing  only  the  first  J  {or  J  4- 
1  if  signal  is  present)  dominant  eigenvectors,  where  7  is  the  number  of  independent  narrowband  Jammers. 
Hence,  they  often  assume  knowledge  of  7  or  require  that  7  be  estimated  by  semi-heuristic  techniques  such 
as  the  AIC  or  MDL  [7].  However,  these  techniques  are  not  derived  via  fundamental  criteria  such  as  a 
ML  estimate,  as  is  the  SMI.  Also  they  are  often  limited  in  their  applicability  or  require  intuitive  rules  to 
determine  the  dimension  and/or  basis  of  the  subspace.  In  [8,9],  we  presented  a  new  technique  that 
provides  typically  similar  performance  as  the  heuristic  techniques,  and  does  not  require  a  priori 
knowledge  of  7.  The  new  fast  maximum-likelihood  (FML)  technique  assumes  only  knowledge  of  the 
receiver  thermal  noise  level;  convergence  MOE  is  similar  to  many  of  the  fast  converging  heuristic 
techniques,  e.g.  in  a  narrowband  (NB)  Jamming  scenario,  convergence  time  is  on  the  order  of  twice  the 
number  of  NB  Jammers.  The  technique  also  works  for  any  external  interference  environment,  for  example 
for  wideband  Jammers  and  clutter  without  requiring  modification. 

An  adaptive  canceller  is  a  particular  form  of  adaptive  linear  processor  (ALP).  In  general,  the  ALP 
assumes  that  the  desired  input  signal  has  the  form  of  an  W-length  steering  vector,  s,  where  N  is  the 
number  of  sensor  inputs  to  the  ALP  and  it  is  assumed  s's  =  1  ('  denotes  conjugate  transpose).  An 
adaptive  canceller  assumes  that  s  =  (1  0  0  . . .  0)"^  =  ij  (where  T  denotes  transpose)  and  the  adaptive 
weight  on  the  main  channel  equals  one.  For  this  form,  it  is  seen  that  all  of  the  desired  signal  energy  is 
assumed  to  be  in  the  first  input  of  the  ALP.  The  first  input  is  called  the  main  channel  and  the  other  inputs 
are  called  auxiliary  channels.  For  example,  a  common  implementation  of  an  adaptive  canceller  is  the 
adaptive  sidelobe  canceller  (ASLC)  of  a  radar  system.  For  this  configuration  the  main  channel  is  the 
output  of  the  radar’s  main  antenna  and  the  auxiliary  channels  are  formed  from  much  physically  smaller 
auxiliary  antennas  which  are  in  close  proximity  of  the  main  antenna. 
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If  digital  processing  is  employed,  any  ALP  can  be  transformed  into  an  adaptive  canceller 
configuration.  This  is  because  the  multiplication  of  the  A-inputs  to  the  ALP  by  any  N  x  N  nonsingular 
matrix  does  not  change  the  signal-to-noise  (S/N)  performance  measure  of  the  ALP.  It  can  be  shown  that 
there  always  exists  a  nonsingular  matrix  A  such  that  /(s  =  Iq. 

One  of  the  advantages  of  an  adaptive  canceller  configuration  is  that  it  can  be  laid  out  functionally  in 
a  highly  numerically  efficient  parallel/pipelined  cascaded  signal  processing  architecture  as  illustrated  by 
the  generic  cascaded  canceller  (GCC)  pictured  in  Fig.  1 .  For  example,  it  is  known  that  the  SMI  algorithm 
for  s  =  Iq  can  be  replicated  exactly  by  using  the  Gram-Schmidt  (GS)  canceller  configuration  [10,11].  For 
this  configuration,  canceller  inputs  (from  the  right  in  the  figure)  are  sequentially  orthogonalized 
(decorrelated)  with  respect  to  other  canceller  inputs.  Data  is  inputted  as  an  A  x  AT  block  of  data  where 
K  is  the  number  of  samples  per  input  channel  used  to  calculate  the  weights.  After  the  A  —  1  weights  are 
calculated  at  the  1st  level  of  GS  canceller,  the  N  x  K  data  block  is  weighted  properly  and  passed  as  an 
(N  —  1)  X  K  data  block  to  the  second  level  for  processing.  Thereafter  a  new  N  x  K  input  data  block 
can  be  processed  by  the  1st  level.  This  sequential  processing,  whereby  input  data  blocks  are  processed 
and  passed  from  level-to-level,  results  in  an  enhanced  processing  throughput  rate.  The  advantage  of  the 
cascaded  canceller  configuration  is  not  that  it  requires  significantly  less  numerical  operations  (which  could 
be  implemented  as  software)  but  that  it  can  be  laid  out  in  a  highly  parallel/pipelined  structure  in 
hardware. 

In  this  paper,  we  develop  what  we  call  a  fast  cascaded  canceller  (FCC).  The  FCC  will  have  almost 
identical  convergence  performance  of  the  FML  canceller  for  interference  scenarios  where  the  number  of 
jammers  is  upperbounded.  By  this  we  mean  that  if  the  number  of  Jammers  is  less  than  a  certain  number, 
the  convergence  performance  of  the  FCC  and  FML  are  almost  identical.  When  the  number  of  Jammers 
exceeds  this  upperbound,  there  is  a  graceful  degradation  in  FCC  performance  with  respect  to  the  FML 
as  the  number  of  Jammers  increase.  However,  it  was  observed  for  A  <  12  and  representative  Jamming 
scenarios  where  the  number  of  Jammers  varied  from  1  to  A  —  1  (normally  the  upper  limit  on  Jammers 
such  that  the  Jamming  can  be  effectively  eliminated),  that  the  FCC’s  average  convergence  performance 
is  within  1  dB  of  the  FML’s;  i.e.  if  the  FML’s  performance  is  —3  dB  from  the  optimum  for  a  given 
number  of  input  sample  vectors,  then  the  FCC’s  performance  is  no  greater  than  —4  dB  from  optimum. 
For  the  FCC  development  we  assume  (as  was  done  for  the  FML)  that  the  unknown  interference 
covariance  matrix  has  the  structure  of  an  identity  matrix  (associated  with  system  noise)  plus  an  unknown 
psdh  matrix  (associated  with  the  external  noise). 

II.  Review  of  FML  Canceller 

The  ML  estimate  for  the  A  by  A  covariance  matrix  that  results  under  the  Gaussian  assumption  and 
AT  >  A  is  the  sample  covariance  matrix  (SCM).  The  convergence  MOE  is  on  the  order  of  AT  =  2A.  We 
derived  in  [8,9]  the  ML  estimate  assuming  a  structured  covariance  matrix  of  the  form 

A  =  c/^  +  Aq, 

where  Rq  is  positive  semi-definite  Hermitian  (psdh),  K  <  N,  the  rank  of  Aq  =  Mq  <  A,  and  Ifj  is  the  A  x 
A  identity  matrix.  We  assume  that  the  thermal  noise  level  c  is  known  and  that  Mq  is  unknown.  Without 
loss  of  generality,  we  take  c  =  1  in  Eq.  (1).  In  much  of  the  adaptive  literature,  the  thermal  noise  level 
has  been  assumed  unknown.  Justification  for  this  assumption  can  be  traced  to  several  papers  [12,13]. 
Certainly,  the  thermal  noise  level  is  unknown  in  many  time-series  problems  such  as  in  the  analysis  of  sun¬ 
spot  data.  However,  in  systems  operating  at  the  microwave  frequencies,  the  thermal  noise  level  is 
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dominated  not  by  the  unknown  external  thermal  noise  but  by  the  receiver  thermal  noise  [14,  p.  1.4], 
Hence,  the  thermal  noise  power  level  is  known  a  priori  for  many  practical  interference  cancellation 
problems.  Under  these  conditions  we  have  found  a  fast  converging  ML  solution  for  the  case  K  <  N. 

The  ML  solution  for  the  structured  matrix  R  -  ly  +  /?„  has  beeii-derived  in  previous  published  work 
for  the  case  K  >  N  and  the  SMC  is  nonsingular.  In  [8,9J  we  showed  that  a  solution  also  surprisingly 
exists  for  K  <  N.  In  fact,  a  solution  can  be  found  even  for  a  single  sample.  There  are  only  a  few  cases 
[15]  of  constrained  covariance  matrices  where  a  ML  solution  exists  for  such  small  sample  cases,  and  most 
of  these  are  not  relevant  to  existing  problems.  Our  solution  is,  however,  directly  relevant  to  many  systems 
that  operate  in  the  microwave  region,  since  the  thermal  noise  is  dominated  by  the  receiver  noise,  of  which 
the  statistics  can  be  assumed  known. 

Let  R  =  Z*  Z^/ K,  where  Z  is  the  Nhy  K  sample  data  matrix  and  R  is  the  SCM.  The  K  columns 
of  Z  are  assumed  to  be  zero-mean,  independent,  and  identically  distributed  (idd)  A'-length  complex 
Gaussian  random  vectors  with  covariance  matrix  R.  The  individual  elements  of  Z  are  complex  circular 
random  variables  (i.e.  the  real  and  imaginary  parts  are  idd).  The  n’th  (n  =  1,  ■ . N)  row  of  Z,  denoted 
by  zl  (a  AT-length  vector),  represents  the  sampled  data  on  the  n’th  sensor.  The  data  in  a  given  column 
of  Z  is  assumed  to  be  time  (or  range)  coincident.  Hence,  each  Af-length  random  vector  (a  column  of  Z) 
is  often  called  a  snapshot  of  input  data. 

The  joint  probability  density  function  (PDF)  of  the  data  under  the  Gaussian  assumption  and  R 
nonsingular  is 


p{Z\  /?„)  =  exp  -  Tr(KR-'R),  (0 

where  R  =  ly  -\-  Rq,  and  |  "I  and  Tr(-)  denote  determinant  and  trace,  respectively.  The  ML  estimate  (if 
it  exists)  for  the  covariance  matrix  is  given  by 

R^^^  =  arg  min  In  |/?|  +  Tr(7?''.R).  (2) 

R 

It  is  known  that  the  ML  solution  exists  for  K  >  N  and  the  SCM  is  nonsingular  [16],  Consider  an 
eigenvalue  decomposition  (evd)  of  the  SCM,  R  =  <I>A4>',  where  A  is  an  x  //diagonal  matrix  with 
diagonal  entries  X,  >  X,  . . .  >  X^  that  are  eigenvalues  of  R  and  is  an  //  x  //unitary  eigenmatrix. 
Denote  M  as  the  number  of  eigenvalues  that  are  greater  than  one.  Set  A^  =  Diag 
(X,,  Xj^,  1,  1,  ...,  1)  where  Diag  (•)  denotes  a  diagonal  matrix  with  elements  and  ordering 

specified  by  its  arguments.  The  ML  estimate  for  R  exists  for  K  >  //and  is  given  by  [16]. 

R  =  4'A 

^'ML  '*'■‘10^  • 

In  [8,9],  we  showed  that  this  same  estimate  can  be  used  for  I  <  AT  <  //.  This  result  is  stated  as  the 
following  theorem. 

Theorem  1:  Under  the  condition  R  =  +  Rq,  where  Rq  is  a  psdh  matrix,  the  ML  solution, 

R  =  arg  min^  In  |/?|  +  Tr(/?''A),  is  given  by 

R  =  4>A  (3) 

for  any  AC  >  1,  and  $  and  Aq  are  as  defined  above. 
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The  FML  N-length-vector  weight,  denoted  by  ,  is  proportional  to  Rml^*  vvhere  s  is  the  N- 
length  steering  vector  of  the  desired  signal.  For  the  canceller  configuration,  s  =  I,,.  In  addition  tor  the 
canceller  configuration,  the  first  element  of  is  constrained  to  be  equal  to  one.  If  a",!,',  is  the  (1,1) 

element  of  ^mi.  then 


w, 


(4) 


We  will  functionally  denote  the  FML  canceller  as  seen  in  Fig.  2.  Ffere  we  also  show  N  /f-length  input 
vectors,  {n  =  1,  2,  ...,  AO .  The  input  z,  is  designated  as  the  main  channel  and  (n  =  2,  3,  N) 
are  called  the  auxiliary  channels.  The  output  of  the  FML  structure  is  the  ALlength  vector  weight, 


We  apply  the  weight  Wpi^^  to  data  (a  snapshot)  that  is  statistically  independent  of  the  data  that  was 
used  to  calculate  Wp^^  and  denote  this  snapshot  as  z.  This  will  be  referred  to  as  nonconcurrent 
processing.  The  snapshot  of  the  data,  z,  is  also  identically  distributed  as  the  snapshots  of  the  data  that  are 
used  to  calculate  Wp,^,  .  We  apply  vVpj^,p  to  z  to  form  the  scalar  output  residue 


We  will  examine  and  compare  the  statistical  characteristics  of  r.  with  the  FCC  output  residue  in  the 
subsequent  analysis. 


III.  Fast  Cascaded  Canceller 


The  general  functional  structure  of  a  cascaded  canceller  was  shown  in  Fig.  1.  The  fundamental 
building  block  of  the  structure  is  the  2-input  canceller  whereby  data  in  the  right  input  of  the  2-input 
canceller  is  weighted  by  a  complex  scalar  and  thereafter  subtracted  from  the  left  input.  For  example  if 
u  and  V  denote  the  A'-length  vectors  of  the  left  and  right  inputs,  w  denotes  the  complex  scalar  weighting, 
and  r  the  AT-length  output  residue  vector,  then  r  =  u  -  wv .  (Note  from  now  on  for  notational  purposes 
we  write  without  loss  of  generality  r  =  u  +  wv  where  w  =  -w).  We  will  associate  the  ij  2-input 
canceller  with  w.^.;  i.e.  the  i,j  2-input  canceller  can  be  found  at  the  /’th  level  and /th  column  of  the  GCC 
seen  in  Fig.  1.  Assume  that  the  input  to  GCC  are  the  N  A'-length  input  vectors  z^(«  =  I,  2,  K) . 
Define  xf  y  =  1,  2,  . . .,  A^  -  /  +  1  as  the  AT-length  output  vector  of  the  iJ  2-input  canceller  with 


Thus 


A'*i) 


-  (fl 


i  =  1,2,  ...,N 
j  =  1,  2,  ...,  N  -  i  ^  1. 


(7) 


The  weight  w.j  is  a  function  of  the  inputs  of  the  i,j  2-input  canceller.  For  example,  for  the  GS  2-input 
canceller 


vv.. 

V 


-(xwl,LixJ^)/  llxf  IP 


where  (|-|  denotes  the  vector  magnitude. 


(8) 


For  the  FCC,  we  use  the  2-input  modified  FML  (MFMLj)  as  a  building  block  as  illustrated  in  Fig. 
3.  The  MFML2  will  be  described  in  detail  in  the  following  section.  The  2-input  MFML2  is  used  to  find 
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the  scalar  complex  weiglit  for  each  /,./  canceller.  However  the  2-input  FML  algorithm  must  be  modified 
because  the  covariance  matrix  of  the  2  inputs  into  the  ij  2-input  canceller  associated  with  the  internal 
thermal  noise  is  no  longer  a  2  x  2  identity  matrix,  A,  except  on  the  first  level.  Furthermore,  it  can  be 
derived  exactly  because  we  know  how  the  input  noises  on  the  various  channels  have  been  weighted  as 
they  traver.se  through  the  cascaded  canceller  structure  up  to  the  iJ  2-input  canceller.  In  order  to  calculate 
the  2x2  internal  noise  covariance  matrix  at  the  2-input  cancellers,  it  is  necessary  to  know  how  each 

input  channel  (/t  =  1 . AO  is  weighted  at  the  output  of  the  iJ  2-input  canceller.  Let  y.j  represent 

the  output  of  the  iJ  MFML,.  We  define  to  be  weighting  on  the  z„  at  the  i,j  output.  Simply 

stated 

y.j  =  E 


It  can  be  shown  the  h'~^{j,n)  can  be  found  reiteratively  as  follows: 

j  =  1,  2,  N  -  i 

+  w.jh''\N  -  /  +  l,rt),  ^  j  2  N 

Initial  condition:  h^'^(j,n)  =  j,  n  =  1,  A^  where  =  1  ify  =  n,  0  otherwise. 


How  the  /2%',/i) ’s  are  used  will  be  described  in  the  next  section. 

If  we  are  performing  noncurrent  processing  then  we  are  interested  in  the  equivalent  weighting  of  the 
cascaded  canceller  structure  illustrated  in  Fig.  3.  We  desire  to  know  the  equivalent  A^-length  weighting 
vector,  denoted  by  passing  the  iV-length  nonconcurrent  data  vector  z  through  the  FCC.  It  is 

straightforward  to  show  that 

1),  2),  ..., 


with  1 ,  1 )  =  1  • 


IV.  2-Input  Modified  FML 


In  this  section,  we  discuss  in  more  detail  the  implementation  of  the  2-input  MFMLj.  As  it  was 
previously  mentioned  the  2x2  internal  noise  covariance  matrix  of  the  inputs  of  the  iJ  2-input  canceller 
is  no  longer  equal  to  /j-  Denote  this  2x2  covariance  matrix  as  R.j.  We  can  derive  in  simple  fashion  this 
covariance  matrix  with  knowledge  of  h'{j,ri)  and  h'{N  -  i  +  \,  n),  (ji  -  1,  2,  ...,  A/).  Set 
h*'’  =  [/i<'^(/,  1),  ....  /i''^{/.A/)f  •  From  (9)  and  the  fact  that  the  internal  noise  power  on 

Z^n  =  1,  ....  A/)  is  1,  it  follows  that 


4(1,  1)  =  hf'hf, 
4(2,  1)  =  h);4,h4 
4(1,  2)  =  R*  (2,  1), 

Rjji^,  2)  =  hv.,,i . 


(12a) 

(12b) 

(12c) 

(12d) 


Let  R.J  be  the  2  x  2  SCM  of  the  iJ  2-input  canceller.  The  elements  of  4  ^re  found  as 
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4(1,  1)  =  xf  xj'V/r, 
4(2,  1)  = 

4(1,  2)  =  R*  {2,  1), 

4(2,  2)  = 


(1 3a) 

(1 3b) 
(13c) 

(13d) 


Let  R.j  denote  the  true  2x2  covariance  of  the  i,j  canceller.  Now  R.j  =  R..  +  R^..  where  R^..  is  some 
2x2  psdh  matrix.  In  order  to  use  the  FML  methodology,  the  internal  noise  covariance  matrix  must  be 
the  identity  matrix.  Thus  we  whiten  the  internal  noise  on  the  input  to  the  ij  2-input  canceller.  Let 
4  “  Cholesky  decomposition  of  R.j  where  A.j  is  a_2  x  2  upper  triangular  matrix  (for  a 

2x2  hermitian  matrix,  this  is  easily  found).  We  form  a  SCM,  R.j  as 

=  C4) 

Let  R.J  be  the  true  2x2  covariance  matrix  associated  with  R.j .  The  2x2  internal  noise  covariance 
matrix  component  of  4  is  4-  ^^^w  estimate  R'.j  via  its  ML  estimate  (see  Theorem  1).  Call  this 

estimate  4  ■  The  ML  estimate  of  4  is  then 


<  =  4  4' 


(15) 


The  w.j  is  found  by  using  (4).  It  can  be  shown  that 


w.. 

u 


In  order  to  find  it  may  be  necessary  to  perform  the  evd  of  R.j.  It  is  pointed  out  that  for  a  2  x 
2  matrix  this  is  very  straightforward  and  that  one  can  derive  closed  form  solutions  for  the  two 
eigenvalues/eigenvectors.  In  fact,  because  the  sum  of  the  outer  products  of  the  2  eigenvectors  is  equal 
to  In,  it  will  be  found  that  only  one  of  the  eigenvectors  must  be  found.  _Furthermore,  finding  this 
eigenvector  is  not  always  necessary.  Let  \  and  X2  be  the  two  eigenvalues  of  /?..  with  \  >  If  > 
1  or  Xj  <  1  then  it  is  not  necessary  to  compute  the  eigenvector.  This  is  because  if  X2  > 
and  if  X.  < 


<(2.  1) 

2) 


(16) 


1  then  4" 


1  then^""  =  R.J 


R.J.  It  was  noted  in  our  simulations,  assuming  correlated  interference  was 
present  at  the  input  of  the  FCC,  that  the  2-input  cancellers  on  the  upper  levels  of  the  FCC  implement  the 
weights  using  Rl^  =  R.j  (which  is  the  GS  weight  implementation).  After  the  interference  had  been 
significantly  reduced  (somewhere  in  the  middle  levels  of  the  FCC),  the  weights  were  calculated  using  the 
evd.  At  the  bottom  level  of  the  FCC  the  weights  were  calculated  using  =  R.j. 


VI.  Results  and  Discussion 


In  this  section,  we  present  representative  simulation  results  that  compare  the  FCC  and  FML 
convergence  performances.  We  shall  see  that  for  interference  scenarios  where  the  number  of  jammers 
is  upperbounded  that  the  FCC  and  FML  convergence  performances  are  almost  identical.  We  also  compare 
the  FCC  and  FML  convergence  performance  with  a  canceller  that  we  call  the  pseudo  sample  matrix 
inversion  (PSMI)  canceller.  For  K  >  N,  the  PSMI  is  exactly  the  SMI  algorithm.  However  for  K  <  N, 
since  the  SMI  algorithm  is  not  applicable  (the  SCM  is  singular),  we  use  an  alternate  procedure  which  we 
now  describe  to  find  the  canceller  weighting  vector. 
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Let  Z.,  be  the  (N  -  \)  x  K  data  matrix  associated  with  auxiliary  channels  fc,  . . z^)  and  be  the 
AT-lensith  column  vector  associated  with  the  main  channel  (z,).  We  can  find  a  weighting  vector  of  length 
N  —  1,  such  that 

(17) 


^7  r 

2a  = 


In  fact  for 
pseudoinverse: 


—  1,  there  are  an  infinite  number  of  solutions  for  w^.  We  solve  the  w^,  using  the 


vv,  =z;(z;z;)-'z„ 


(18) 


The  pseudoinverse  solution  minimizes  ||  ||  with  respect  to  all  possible  solutions  of  (17).  The  AAlength 

weight  vector  (1,  -w^)^  is  the  solution  that  we  use  for  the  PSMI  when  K  <  N  —  1.  For  K  >  N,  the 
SMI  algorithm  is  used  to  find  w^.  As  with  the  FML  and  FCC  the  PSMI  is  implemented  using 
nonconcurrent  data. 

Although  there  are  many  configurations  that  can  be  considered  for  interference  cancellation  problems, 
for  simplicity,  we  assume  that  an  N  element  array  of  identical  antenna  elements  exists  such  that  the 
jammer  vectors  have  the  following  form 

(1,  exp  (j^p,  exp  (2j0p,  ...,  exp  {(N  -  l)J^p)^, 

where  j  =  /^T.  The  main  channel  is  the  left-most  array  element  and  the  remaining  N  -  1  elements  are 
the  auxiliary  channels.  We  model  the  external  noise  environment  via  its  input  covariance  matrix.  The  N 
X  N  covariance  matrix  associated  with  J  narrowband  Jammers  can  be  represented  as  =  (r  where 


=  E  [j('^ " 


(20) 


<  =  l 


</)  (/  =  1,2,  . . .,  J)  is  the  phase  angle  associated  with  the  I’th  jammer,  and  a~  is  the  jammer  power  of 
the  z’th  jammer  normalized  by  the  internal  noise  level.  The  5,„„  contribution  to  the  covariance  matrix  is 
associated  with  the  internal  thermal  noise  power. 

In  order  to  compare  different  techniques,  we  will  plot  the  normalized  average  signal-to-interference 
ratio  (SIR),  p,  which  is  defined  as  the  averaged  output  SIR  ratio  for  a  given  technique  divided  by  the 
optimal  SIR  that  can  be  achieved  when  the  optimal  weight  is  used.  Interference  as  defined  here  includes 
all  unwanted  interference  including  thermal  noise.  The  optimal  SIR  is  known  to  be  SIR^  =  IqR  '1 
Averaged  over  the  Monte  Carlos,  p  is  defined  by 

|wIl„P 


0  ■ 


1 


MC 


m  =  l 


w'Rw  SIR 

m  in  opi 


(21) 


where  w„  is  the  random  weight  vector  associated  with  a  given  canceller  technique  (FML,  FCC,  or  PSMI) 
indexed  by  the  Monte  Carlo  number  and  MC  denotes  the  number  of  Monte  Carlos.  The  kernel  of  the  sum 
seen  in  (21)  represents  the  normalized  instantaneous  SIR  for  a  given  Monte  Carlo.  The  expected  value 
of  the  nonconcurrent  data  has  been  taken  and  is  exemplified  by  the  R  term  in  the  denominator  of  the 
kernel. 

For  some  of  our  simulation  results,  we  fix  the  jammer  angles  and  powers  and  generate  MC 
realizations  of  input  data.  For  the  rest  of  our  simulations,  we  make  the  jammer  angles  and  powers  random 
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variables  and  generate  MC  realizations.  For  this  latter  case,  S1R„,„  is  a  function  of  the  jammer  angles  and 
powers  for  a  given  realization.  Hence  in  (21),  SIR,,,,  will  vary  as  m,  the  Monte  Carlo  index. 

We  now  show  a  number  of  simulation  results  for  representative  interference  scenarios.  In  all  cases 
the  number  of  input  channels  to  the  canceller,  N,  is  20  and  the  number  of  Monte  Carlos  (MC)  equals 
100.  In  Figs.  4-8,  canceller  performance,  normalized  average  S/I,  is  plotted  vs.  the  number  of 
independent  snapshots  K  for  the  canceller  configurations:  FCC,  FML,  and  PSMI  and  various  interference 
scenarios.  The  maximum  K  was  chosen  to  be  equal  to  40  (=2N)  which  is  normally  the  K  chosen  for 
“good”  convergence  performance  (-3  dB  on  average  from  the  optimal)  for  the  SMI.  We  see  from  Figs. 
1-4  that  for  7  <  4,  the  FCC  and  FML  convergence  performances  are  almost  identical.  The  PSMI 
canceller  performance  is  notably  inferior  to  the  FCC  and  FML  over  much  of  the  range  of  K.  It  is 
interesting  to  note,  however,  that  for  I  <  K  <  N,  a  maximum  occurs  in  the  PSMI  convergence 
performance,  which  in  a  number  of  cases  is  not  much  less  than  the  FCC  and  FML  performance. 

We  see  from  Fig.  8  that  for  /  =  5  the  FCC  convergence  performance  is  noticeably  less  than  FML’s. 
It  was  observed  over  a  number  of  cases  (not  presented  here)  that  the  FCC’s  convergence  performance 
for  a  given  Jamming  scenario  (for  certain  Jammer  angles)  was  close  to  the  FML’s  up  to  a  certain  number 
of  Jammers.  After  the  required  number  of  Jammers  had  been  reached,  the  convergence  performance 
decreased  monotonically  with  the  increase  in  the  number  of  Jammers.  However,  it  was  observed  (these 
plots  are  not  presented  here),  for  A/  <  12  and  representative  Jamming  scenarios  where  the  number  of 
Jammers  varied  from  1  to  A  —  I  (normally  the  upper  limit  on  Jammers  such  that  the  Jamming  can  be 
effectively  eliminated),  that  the  FCC’s  average  convergence  performance  is  within  1  dB  of  the  FML’s; 
i.e.  if  the  FML’s  performance  is  -3  dB  from  the  optimum  for  a  given  number  of  input  sample  vectors, 
then  the  FCC’s  performance  is  no  greater  than  —4  dB  from  optimum. 

In  Figs.  9-12,  we  plot  the  normalized  average  S/I  averaged  over  a  number  of  Jammer  angles  and 
power  realizations  vs.  K.  For  a  given  number  Jammers,  7,  for  each  Monte  Carlo  we  choose  7  Jammer 
angles  independently  which  are  uniformly  distributed  on  [0°,  360°].  In  addition,  we  choose  each 
Jammer’s  power  from  a  uniform  distribution  on  [15  dB,  40  dB].  From  Figs.  9-12  we  observe  that  the 
FCC’s  performance  is  almost  exactly  the  same  as  the  FML’s  for  7  <  5.  For  7  >  5,  there  is  a  gradual 
decline  in  performance. 

The  sampled  standard  deviation  (s.d.)  for  a  few  of  the  cases  presented  in  Figs.  4-12  are  shown  in 
Figs.  13-16.  Here  we  observe  that  the  FCC’s  s.d.  is  similar  to  the  FML’s  when  the  FCC  convergence 
performance  is  close  to  the  FML’s  and  that  the  FCC’s  and  FML’s  s.d.  is  moderately  better  than  the 
PSMI’s  for  K  approaching  2N  (or  40).  In  fact  one  might  argue  that  one  reason  for  using  the  FML  (or  the 
FCC  when  its  convergence  performance  is  almost  identical  to  the  FML)  rather  than  the  SMI  for  K  —  2N 
is  that  the  FCC  or  FML  output  s.d.  of  the  output  residue  is  noticeably  better  than  the  SMI. 

The  number  of  floating  point  multiplication  operations  (FPMOP)  associated  with  finding  the  canceller 
weight  for  the  FCC  is  approximately  .5N^{3K  -f  20).  For  the  FML  canceller  implemented  with  full 
singular  value  decomposition  (svd)  (all  singular  values  and  eigenvectors  calculated  via  svd),  the  number 
of  FPMOPs  is  approximately  3N^K  +  4NK^ .  For  the  FML  canceller  implemented  with  partial  svd  (only 
the  eigenvalues  greater  than  1  and  their  associated  eigenvectors  need  be  calculated)  the  number  of 
FPMOPs  =  O(KNJ)  (this  assumes  a  NB  Jamming  scenario).  For  the  SMI  canceller, 
FPMOP  *=  .33N^  +  KN-.  Hence  the  number  of  FPMOP’s  needed  to  implement  the  FCC  in  software 
is  roughly  equivalent  to  either  the  FML  or  SMI  implementations.  However,  as  it  was  pointed  out  in  the 
introduction,  the  advantage  of  the  FCC  is  not  its  software  numerical  efficiency  but  its  hardware  numerical 
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efficiency.  The  algorithm  can  be  laid  out  functionally  using  a  highly  parallel/pipeline  architecture.  Thi.s 
structure  is  ideal  for  efficiently  block  proces.sing  input  data  blocks  that  are  sequentially  updated  at  each 
time  step  (or  some  other  dimension). 

VII.  Suniniary 

A  fast-converging,  highly  parallel/pipeline  cascaded  canceller  (FCC)  has  been  developed  which  has 
convergence  performance  almost  identical  to  the  fast  maximum  likelihood  (FML)  canceller  [8,9]  for 
restricted  jamming  scenarios.  For  narrowband  jamming  scenarios,  it  has  been  shown  that  the  FCC 
convergence  performance  is  similar  to  the  FML’s  when  the  number  of  jammers  is  below  some 
upperbound.  However,  it  was  observed  for  /V  <  12  and  representative  jamming  scenarios  where  the 
number  of  jammers  varied  from  1  to  V  -  1  (normally  the  upper  limit  on  jammers  such  that  the  jamming 
can  be  effectively  eliminated),  that  the  FCC’s  average  convergence  performance  is  within  1  dB  of  the 
FML’s;  i.e.  if  the  FML’s  performance  is  -3  dB  from  the  optimum  for  a  given  number  of  input  sample 
vectors',  then  the  FCC’s  performance  is  no  greater  than  -4  dB  from  optimum.  For  both  the  FML  and 
FCC  developments  it  is  assumed  that  the  unknown  interference  covariance  matrix  has  the  structure  of  an 
identity  matrix  plus  an  unknown  positive  semi-definite  Hermitian  (psdh)  matrix.  The  identity  matrix 
component  is  associated  with  the  known  covariance  matrix  of  the  system  noise  and  the  unknown  psdh 
matrix  is  associated  with  the  external  noise  environment.  For  narrowband  jamming  scenarios  with  J 
jammers  where  J  is  less  than  some  upperbound,  it  was  shown  via  simulation  and  analysis  that  the  FCC 
and  FML  converge  on  the  average  —  3  dB  below  the  optimum  in  about  2J  independent  sample  vectors 
per  sensor  input.  Both  the  FCC  and  FML  converged  much  faster  than  the  SMI  algorithm. 
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Fig.  3  —  Fast  cascaded  canceller. 
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Normalized  S/I  vs.  no.  of  independent  snapshots,  K\  N  -  20,  J  2,  dj  40  ,  60  ,  P 
=  100,  three  canceller  configurations:  fee,  fml,  psmi. 


normalized  average  S/I 


no.  of  independent  snapshots,  K 

Fig.  6  —  Normalized  S/I  vs.  no.  of  independent  snapshots,  K-,  N  =  2Q,  J  -  3,  dj  —  20°,  40°,  60°,  Pj 
=  30  dB,  MC  =100,  three  canceller  configurations:  fee,  fml,  psmi. 
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normalized  average  S/I 


normalized  average  S/I 


jammer  angles  invariant  over  monte  carlos 


Fig.  8  —  Normalized  S/I  vs.  no.  of  independent  snapshots,  K;  N  =  20,  J  =  5,  dj  =  20°,  30°,  40°,  50°, 
60°,  Pj  =  30  dB,  MC  =  100,  three  canceller  configurations;  fee,  ftnl,  psmi. 


18 


Fig.  10  —  Normalized  S/I  vs.  no.  of  independent  s 
distributed  [0,  360°],  jammer  powers  uniformly 
configurations:  fee,  fml,  psmi. 


normalized  average  S/! 


averaged  ove 


Fig.  11  —  Normalized  S/I  vs.  no.  of  independent 
distributed  [0,  360°],  jammer  powers  uniformlj 
configurations:  fee,  fml,  psmi. 
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Fig.  12  —  Normalized  S/I  vs.  no.  of  independent  snapshots,  K]  N  =  20,  J  =  10,  jammer  angles 
uniformly  distributed  [0,  360°],  jammer  powers  uniformly  distributed  [15,  40]  dB,  MC  =  100,  three 
canceller  configurations:  fee,  fml,  psmi. 


sampled  standard  deviation 


jammer  angles  invariant  over  monte  carlos 


no.  of  independent  snapshots ,  K 


Fig.  14  —  Sampled  standard  deviation  vs.  no.  of  independent  snapshots,  K;  N  —  20,  J  —  5,  6j 
30°,  40°,  50°,  60°,  Pj  =  30  dB,  MC  =  100,  three  canceller  configurations:  fee,  ftnl,  psmi. 


Fig.  16  —  Sampled  standard  deviation  vs.  no.  of  independent  snapshots,  K-,  N  =  20,  J  =  7,  jammer 
angles  uniformly  distributed  [0,  360°],  jammer  powers  uniformly  distributed  [15,  40]  dB,  MC  =  100, 
three  canceller  configurations;  fee,  fml,  psmi. 


